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We develop practical schemes for the measurement of the covariance matrix for intrinsic angular- 
momentum variables in quantum optics. We particularize this approach to two-beam polarimetry 
and interferometry, as well as to ensembles of two-level atoms interacting with classical fields. We 
show the practical advantages of noisy simultaneous measurements. 



I. INTRODUCTION 

Angular-momentum variables represent basic observ- 
ablcs both in classical and quantum optics, especially 
in three fundamental areas: polarization, interferometry, 
and light-matter interaction [J, 0, H H B, @, 0, B S El • 
For definiteness, throughout we focus on intrinsic (not 
orbital) angular momenta. This is for example the case 
of the Stokes parameters, which provide a complete ac- 
count of second-order (in complex amplitudes) statisti- 
cal properties of two-mode polarization and interference. 
Moreover, spin operators are basic in atomic physics such 
as in the case of ensembles of two-level atoms described 
individually as spin 1/2 systems. 

The second-order statistics of angular-momentum vari- 
ables are crucial in diverse areas. This is the case with 
quantum metrology, where angular-momentum statistics 
determine the ultimate limit to the resolution of in- 
terferometric and spectroscopic measurements 0, d, 
Moreover, angular-momentum covariance matrices en- 
ter in the analysis of many-body entanglement in 
continuous- variable polarization entanglement 0], and 
for light-mediated detection of atomic-spin correlations 
1- 

Recently we have proposed an SU(2)-invariant char- 
acterization of angular-momentum fluctuations via the 
diagonalization of the covariance matrix Invariance 
under SU(2) transformations is a desirable property since 
two states connected by a deterministic SU(2) transfor- 
mation should be statistically equivalent. Similar invari- 
ance ideas are at the heart of current investigations about 
the coherence between classical vectorial waves [l2| ■ 

In this work we develop simple practical schemes to 
determine experimentally the angular-momentum covari- 
ance matrix of a given system in an unknown state. We 



particularize the method to diverse optical two-mode po- 
larimetric and interferometric configurations, as well as 
to ensembles of two-level atoms. It is worth stressing that 
this analysis applies equally well to quantum and classi- 
cal optics. In the classical domain the situation is much 
more simple since in principle one can always perform as 
many simultaneous measurements as desired of any set of 
observables in accurate copies of the original beam pro- 
vided by beam splitting, for instance. This idea can be 
fruitfully translated to the quantum domain in the form 
of noisy simultaneous measurements of noncommuting 
angular-momentum components. 

In Sec. II we recall basic definitions and results. In 
Sec. Ill we present a basic scheme for the measure- 
ment of angular-momentum covariance matrices which is 
particularized to polarimetric, interferometric, and spec- 
troscopic situations. In Sec. IV we present an inter- 
ferometric noisy simultaneous measurement of angular- 
momentum components providing a simple and exact 
practical determination of the covariance matrix with a 
single experimental configuration. In Sec. V we consider 
the bright limit in which the angular-momentum covari- 
ance matrix becomes a quadrature (or position- linear mo- 
mentum) covariance matrix. 



II. DEFINITIONS 

A. Definition and two-mode realization 

Let us consider arbitrary dimcnsionless angular mo- 
mentum operators j* = (ji, J2,ja), where the superscript 
t denotes matrix transposition. In quantum physics they 
are defined by the fulfillment of the commutation rela- 
tions 



'Electronic address: A.Rivas@herts.ac.uk 

t Electronic address: 

URL: http://www.ucm.es/info/gioq 



alluis@fis.ucm.es 



\jk,jt] 



UoJ] 



(2.1) 



2 



where tk,e,n is the fully antisymmetric tensor with £1.2, 3 
1 and j is defined by the relation 



jo (jo + 1) 



(2.2) 



For the sake of completeness we take into account that 
jo may be an operator. This is the case of two-mode 
realizations where jo is proportional to the number of 
photons. In classical optics the situation is similar by 
replacing commutators by Poisson brackets and Eq. (|2.2j) 
by P = jo (with (j) 2 < (jo) 2 , where the brackets denote 
an ensemble average). 

In quantum and classical optics two-mode realizations 
of angular momentum play a relevant position in po- 
larimetry and interferometry. In the quantum case, de- 
noting by ai ; 2 the complex amplitudes operators of two 
field modes, we get that 
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J2 = 



a\a 2 



J3 = 



a\a\ 



a\a2 



satisfy Eqs. (|2.1[) and (|2.2p [13| . where the superscript f 
denotes Hermitian conjugation. In the classical domain 
0-1.2 ar e classical amplitudes so that Hermitian conjuga- 
tion a\ 2 is replaced by complex conjugation a\ 2 . 

In polarimetry these are essentially the Stokes vari- 
ables. The normalized vector (j) / (jo) defines the 
Poincare sphere as a suitable representation of polariza- 
tion states and transformations 0, [l(|. These are also 
basic variables in two-beam interferometry. For example, 
the third component j3 is proportional to the difference 
of the number of photons between two modes, while ji,i 
express the coherence between the interfering beams. 



B. Covariance matrix 

The complete second-order statistics of j is contained 
in the 3x3 real symmetric covariance matrix M 

1 



M k,i = 77 ((jkje) + (jejk)) - (jk) (je) , 



(2.4) 



with M* = M and M* = M. The alternative definition 
M'f. £ = (jkje) — (jk) (je) is identical to M in the classical 
case, while in the quantum domain it provides a com- 
plex Hermitian matrix that contains essentially the same 
information as M (Tlj . 

The covariance matrix M allows us to compute the 
variance (Aj u ) 2 of an arbitrary angular- momentum com- 
ponent j u — u ■ j, where u is any unit real vector, 

:J u'Mtt, 



(Aj.) 



(2.5) 



as well as the symmetric correlation of two arbitrary com- 
ponents j u — u - j and j v — v ■ j, where u and v are unit 
real vectors, 
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((jujv) + (jvJu)) - (ju){jv) = v f Mu = u l Mv. (2.6) 



Since M is real and symmetric, the transformation that 
renders M diagonal is a rotation matrix Rd 



M = Rl 



(AJi 







(AJ2) 2 
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(AJ3) 2 . 



R d . (2.7) 



The eigenvalues of M, (AJ fc ) 2 , k = 1,2,3, are the vari- 
ances of the components Jk — Uk ■ j , where Uk are the 
three real orthonormal eigenvectors of M 



Mu k = {AJ k ) 2 u k . 



(2.8) 



Following standard nomenclature in statistics we refer to 
J and AJ as principal components and principal vari- 
ances, respectively. We stress that both J and AJ de- 
pend on the system state. The principal variances pro- 
vide an SU(2) invariant characterization of angular mo- 
mentum fluctuations II ill . 



C. SU(2) invariance 

Throughout, by SU(2) invariance we mean the statisti- 
cal equivalence between states connected by unitary de- 
terministic transformation generated by j 



U — exp (i9u ■ j) 



(2.9) 



where 9 is a real parameter and it is a unit three- 
dimensional real vector. It can be seen (for example, 
by using the jo = 1/2 representation) that the action of 
U on j is a rotation R of angle and axis u [l4| 



U*jU = Rj, tfj U 



where the 3x3 real matrix R is 



: Jo, 



(2.10) 



Rk,e = {virfcTkU) , (2.11) 



with 



1 3 

U = exp (tW) , V = - 



Uk&k, 



(2-12) 



fc=i 



where k, I = 1,2,3, 01,2,3 are the Pauli matrices and it 
holds that R* = Rr 1 , and W =U~ 1 . 

The SU(2) invariance of principal variances holds be- 
cause under any SU(2) transformation M transforms as 
M — * RMR 1 . Therefore, the covariance matrix RMR 1 
associated with the transformed state has the same prin- 
cipal variances as the covariance matrix M associated 
with the original state. 

In other words, the SU(2) invariance is just the math- 
ematical statement corresponding to the fact that the 
conclusions which one could draw from an angular mo- 
mentum measurement must be independent of which set 
of three orthogonal angular momentum components one 
chooses. 
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In the case of the two-mode bosonic realizations (I2.3|) 
we have 



b = U^aU =Ua, U = exp (i9a^Va) , 



(2.13) 



where U and V are in Eq. (|2.12| and a* = (01,02) and 
b* = (61,62) are the original and transformed complex 
amplitudes, respectively. In this case SU(2) transforma- 
tions describe basic lossless polarization and interference 
elements, such as beam splitters, phase glates^ two^beam 
interferometers, Faraday rotators, etc. 



plates, two- Deai 
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III. PRACTICAL DETERMINATION OF THE 
COVARIANCE MATRIX 

The complete determination of the covariance matrix 
in a given basis of components j* = (ji,j2, J3) can be 
achieved by measurement of the variances of the six op- 
erators 



h±2 = Ui±h) ■ 
ii±3 = ^ (ii ± h) ■ 
ia±3 = ^5 (j'2 ± 33) ■ 



(3.1) 



More specifically, since (respecting the quantum lack of 
commutation) 

il±t = \ 01 + 3e ± Jkjt ± jtjk) , (3.2) 

we get that the nondiagonal matrix elements Mk t e, k £, 
are given by 



M M = ~ [(Aj k+e f - (A Jfc _ £ ) 2 
For the diagonal elements we have 

\2 1 



(3.3) 



3,3 
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(Aji +3 ) 2 + (Aji_ 3 ) 2 + (Aj2 +3 ) 2 



+ (AJ2-3) 2 - (Aji +2 ) 2 - (Aji_ 2 ) 2 , (3.4) 



and similarly for Aji and Aj'2 by cyclic permutations of 
the indices. 

Note that M has just six independent components be- 
cause of reality and symmetry, which agrees with the 
above number of independent measured variances. The 
six operators (|3.ip are not, strictly speaking, independent 
since we have, for example, 



J2+3 — Jl+2 — Jl-3, J2-3 — ]l+2 — Jl+3- 



(3.5) 



Nevertheless, measurement of the components j2±3 is 
necessary to derive all jkji correlations exclusively in 
terms of variances. 

When one of the components of j, say j'3, is a principal 
component, the process is much more simple since we 
know in advance that all the correlations between j'3 and 



FIG. 1: Poincare sphere illustrating polarization states (cir- 
cular at the poles and linear at the equator) and the action 
of the transformations (13.811. 



iij 32 vanish. Then, only four variances are necessary: 
namely Aji, Aj2, Aj'3, and Aji + 2 for example. 

The measured components j p ± q can be related with 
the original ones je by simple SU(2) transformations of 
the form 



±m 



exp(z6> ±m j fc ) . 



(3.6) 



with 8± m = ±7r/m and m = 2, 4, which produce rota- 
tions of angles ir/2 and 7r/4 around the axis j\. This is 
useful because Uk,±4 transform the measurement of the 
components je in the transformed state into the measure- 
ment of the operators j P ± q in the original state, while 
Uk,±2 transform the components je among themselves. 
The proper use of these transformations is illustrated by 
the following particular cases. 



A. Polar imetry 

Polarization states and transformations can be prop- 
erly represented in the Poincare sphere, as illustrated in 
Fig. 1. For defmiteness we consider 01,2 in Eq. (|2.3[) as 
the amplitudes of circularly polarized modes, 
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V2 
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V2 



(3.7) 



where the complex amplitudes of modes linearly 

polarized along the Cartesian axes x and y. As custom- 
ary, the south and north poles in axis j'3 of the Poincare 
sphere represent circularly polarized light, while linear 
polarizations of different azimuths are distributed along 
the equator, with linear polarization along the Cartesian 
axes x and y located at the antipodal points of the axis 
ji (i.e., ji = ±j ). 

The most simple polarization measurement is the mea- 
surement of ji — (a\,a x — a) y ay)/2 as the difference be- 
tween the field intensities after a polarizing beam splitter, 
as illustrated in Fig. 2. In this case the transformations 
(|3.6|) correspond to phase plates and Faraday rotations 
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FIG. 3: Illustration of the interferometric realization of the 
transformations (|3.9[l and (|3. 10|) using symmetric beam split- 
ters (SBS) and phase-difference shifts (PDS) exclusively. 



FIG. 2: Illustration of the scheme for the polarimetric mea- 
surement of the covariance matrix. 

placed before the polarizing beam splitter that transform 
the measurement of j± in the output fields into the mea- 
surement of jk±e in the input fields. More specifically, 

jl±2 = & r ], ±4 ilf / 3,±4, 
il±3 = t4= F4 jl C/ 2,=F4, 

j2± 3 = r: ,r,.. 2 . (3.8) 

The transformations £^,±4 are Faraday rotations produc- 
ing a phase different shift of ±7r/4 between dextro and 
levo circularly polarized modes. This produces a rotation 
of angle ±7r/8 of the azimuth of linearly polarized light, 
which is a rotation of the Poincare sphere of angle ±7r/4 
along the north-south axis. On the other hand, the trans- 
formations C/2,±2 and J72,±4 can be implemented by phase 
plates introducing phase-difference shifts of ±7r/2 and 
±7r/4, respectively, the phase-plate axes forming ±7r/4 
with the Cartesian axes x and y. 

B. Two-beam interferometry 



In Fig. 3 we illustrate how these transformations may 
be implemented with very simple elements such as sym- 
metric beam splitters (SBS) and phase-difference shifts 
(PDS), described by the unitary matrices 

«PDS - (T exp( °_ #) ). (3.1D 

For U2.± m we have the parameters tp = 7r/2, <j> = 9± m /2— 
7r/2, and S — 0, while for Ui.± m the parameters are tp = 
S = 7r/4 and the same <j> = 9± m /2 — tt/2. 



C. Two- level atoms 

In this case the physical situation corresponds to a col- 
lection of N two-level atoms (with ground \g) and ex- 
cited |e) levels) interacting with a classical field E = 
-Eo cos(wi). By assuming that the coupling between 
atoms can be neglected the total Hamiltonian is given 
by the sum of individual Hamiltonians (in units % = 1 
for simplicity), 



Two-beam interferometry can be embedded in this 
same framework by considering that the complex ampli- 
tudes 01,2 represent two interfering modes with the same 
polarization state and propagating along different direc- 
tions. In this case, the simplest measurement is ja, since 
it represents the difference of intensities between the two 
waves 01,2 • Otherwise, the same relations ()3.8[) hold sim- 
ply by the cyclic permutation (1, 2, 3) — * (3, 1, 2) for the 
indices k and i in j k±e , j k , and U k ,± m - 

The transformations (13. 6| represent in general lossless 
beam splitters and phase shifts. In terms of input-output 
relations (|2.13p we get for Ui,±m the following unitary 
matrices relating input and output complex amplitudes 



N 



u 



_ ( cos(0±/2) sin(0±/2) 
1 m ~ I -sin(0±/2) cos(0±/2) 



while for U±± m , 



_ ( cos(0±/2) isin(0±/2) 
" 1 - ±m ~ Usin(0±/2) cos(0±/2) 



(3.9) 



(3.10) 



k=l 



being 



i w o Ik) n 
h k = —ak - U 



cos(ajt), 



where, for each atom k, 

a 3 = \e){e\-\g){g\, a. = a{ = \g){e\ 



(3.12) 



(3.13) 



(3.14) 



are the corresponding Pauli and ladder matrices with 
cr ± = (ai ± «cr 2 )/2. The first term in Eq. (|3~lU)) is the 
free-evolution Hamiltonian for each atom and the second 
one is the coupling with the classical field at dipolar ap- 
proximation. The Rabi frequency Q = (g\d\e) ■ Eq has 
been assumed to be real. On the regimen ui ~ ujq 3> f2 is 
usually used to consider the rotating wave approximation 
by neglecting the counter-rotating terms er + exp(iw£) and 
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<7_ exp(— icot) in the above Hamiltonians 
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LOO _.(*) 

T °3 



2~ 
n 
2" 



exp(iu>t) + <T y +' exp(— iut) 



(ft) 



cr cos(wi) + cr^ sin(wi) .(3.15) 



|V(t)) - |#)) = WW), 



(3.17) 



with \ip(0)) = |^(Q))- The time-evolution equation in 
this picture is 

ij t m)) = [(wo - w)j 3 - fiji] (3.18) 

so the operator U(t2,ti) performing the time evolution 
between t\ = and ti = t in the Schrodinger picture is 

U(t,0) = exp(-iujtj 3 ) exp{-i [(oj - u>)j 3 - Vtjx] t}. 

This operator reduces to a very simple product of SU(2) 
transformations when the radiation field is in resonance 
lu = ujq: 

U Ies (t, 0) = exp(-icj tj 3 ) exp(i^ji). (3.20) 

On the other hand, if the external field is switched off 
(Cl = 0), the free evolution is given by 



Uf ree (t) = exp(-iaj tj 3 ). 



(3.21) 



Therefore, in this case the transformations Uk± m are ob- 
tained by combining time intervals of resonance pulses 
and free evolution, as is used in Ramsey spectroscopy 

As in the interferometric case above, the simplest mea- 
surement is J3 again, this is the difference of populations 
between the two levels of the atoms (nevertheless, see 
Ref. [|| for other light-mediated atomic-spin measuring 
schemes). More explicitly U2,± m can be achieved as 

U2,±m = Uf- [cc (tn/2 — t± m )U Ies (t± m ,0)Uf Iee (t- n /2), 

(3.22) 

with woi±7r/2 — ±tt/2 mod(27r) and, in order to deal 
always with positive time intervals, Slt m = n/m and 
Qt-m — (2m — l)n/m. Similarly Ui_± m can be achieved 
as 

fl,±m = Uf Iee (t2-n- — t± m )U rcs (t± m , 0), (3.23) 
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From Eqs. (|3.12p and (|3.15[) the total Hamiltonian can 
be written in terms of the total angular momentum j = 
£.V fe) /2as 

H = Uoj 3 -tt [ji cos(wi) + j 2 sm(ut)] = u\ [aj j 3 - ilji] U 3 , 

where U 3 = exp(iojtj 3 ) is a rotation around the j 3 axis. 

It is customary to change the picture by the unitary 
transform U 3 = exp(icutj 3 ) in order to remove the time 
dependence of the Hamiltonian 
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FIG. 4: Illustration of the 12-port scheme. 



with LO<st2-K = 27r mod(27r) and the same t± m above. We 
stress that the mod(27r) freedom should be used to obtain 
always positive time intervals. 

In these equations the following relation is useful 

exp(iipj 2 ) = exp exp(iipji) exp (^jsj , 

(3.24) 

which can be derived from the relations in Sec. IIC for 
9 = tt/2 and V — cr 3 /2, so that the rotation matrix in 
Eqs. (|2~TU]) and (j2~TTj) becomes 




(3.25) 



IV. SIMULTANEOUS MEASUREMENTS 

In this section we show that an interferometric noisy 
simultaneous measurement of the three components of j 
provides a simple and exact determination of the covari- 
ance matrix with a single apparatus (similar schemes may 
be developed for the other contexts). The subject of si- 
multaneous measurements of noncommuting observables 
has a long history, being closely related to basic issues of 
quantum theory such as generalized measurements, state 
reconstruction, complementarity, uncertainty relations, 
etc. [jq . In our case we consider the 12-port scheme 
illustrated in Fig. 4 [T3, EH • The two input signal modes 
are a% and 02, while the input modes a±o, 020, a' w and 
a' 2 Q are auxiliary modes always in the vacuum state. 

For definitcness and to simplify formulas let us con- 
sider that beam splitters BS1, BS2, BS3 and BS4 are 
50% with real transmission and reflection coefficients and 
a 7r phase change in the upper side reflections. BS5 and 
BS6 are identical with real transmission t and reflection 
r coefficients, with t 7^ r, and a 7r phase shift in the upper 
side reflections. Finally PS is a it/ 2 phase shift. The re- 
lation between the input and output complex amplitudes 
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In the classical domain the vacuum state implies that 
aio = <220 = a'i = a' 20 = 0, so we have the noiseless 
simultaneous measurement of all the components (|2.3p 
via the detection of the six output intensities Ij = a*dj, 
j = 3, . . . , 8, in the form 

1 8 1 

JO = o ZJ 7 J' Jl = ^2 ( /4 _ /3 ) ' 

]2 = ^{h-h), J3 = ^(h-Ij). (4.2) 

With this we can compute the whole covariance matrix 
M by determining the variances and correlations between 
the output intensities Ij. 

In the quantum case the amplitudes of the auxiliary 
modes et 10 , a 20 , a' 10 and a' 2Q cannot be taken as zero since 
the complex amplitude of the vacuum fluctuates. In other 
words, simultaneous exact measurements of noncommut- 
ing operators are forbidden by commutation relations. 
Nevertheless, it is still possible to extract useful and reli- 
able information from simultaneous noisy measurements. 
To this end let us define the commuting measured observ- 
ables 

31 = J2 (a\ai - 4a 3 J , 

32 = F ( a 6 a 6 - a l a 5) - 

33 = (ala s - a\a^ , (4.3) 

as providing a noisy joint measurement of the operators 
(12. 3p . We do not include jo because this measurement 
is actually exact and noiseless because of conservation 
of total photon number between the input and output 
(the auxiliary modes are in an eigenstate of the number 
operator). 

InRef. [13 it was shown that for the mean values and 
variances we have 

<J> = U) (4-4) 

and 

(Ah) 2 = (Ajif + ^Oo), 



(A] 2 f = (A j2 ) 2 + ^±f 0„), 

(Aj 3 ) 2 = (Aj 3 f + ^{j ), (4.5) 

so that the diagonal terms of the covariance matrix 
Aji,2,3 can be determined simply and exactly from 
Aj'1,2.3 and (jo)- Concerning the nondiagonal terms, it 
can be seen that the following exact relations hold for all 

tidk) = (jkJt) = (■ 3t3k ■) = ^{(jkh+jdk)), (4-6) 

where : : denotes normal ordering. To derive this last 
relation it can be taken into account that a\aua\ai = 
a\a\a^at for k 7^ I in order to express jijk in normal 
order. This is useful since this automatically removes 
the operators of the auxiliary modes in the vacuum state, 
leading to (jejk) — ('■ jejk '■)■ Then, the last equality in 
Eq. (|4.6p can be proved by direct computation. 

Therefore we get that the statistics of j allow one to 
determine the exact mean values and the covariance ma- 
trix for j . The variances of j present an excess of fluctu- 
ations caused by the vacuum in the auxiliary modes that 
can be easily subtracted or compensated. This is partic- 
ularly simple for t = ^/2/3, r = ^/l/3, since in such a 
case 

M = M — (j )l, (4.7) 

where M and M are the correlation matrices for the j 
and j operators, respectively, and J is the 3x3 identity 
matrix. 

It is worth stressing the simplicity of this method since 
it provides complete information via the measurement of 
just four observables (j, jo) instead of the six observables 
of the general method in Eq. (|3.1| . Moreover, the four 
observables (j, jo) are measured in a single experimental 
arrangement. 

The above relations (|4.6j) can be regarded as a cor- 
respondence between classical variables (the outputs of 
measuring j) and quantum mechanical operators j. In 
particular, Eq. (|4.6p is actually an angular momentum 
version of the Wigner and Terlesky-Margenau-Hill 
[20l | correspondences between products of classical vari- 
ables and symmetric operator orderings. 

Similar schemes may be developed for the polarimetric 
context. The same interferometric scheme above is valid 
if the signal modes ai j2 are the polarization components 
a XiV . A more polarimetric scheme allowing the noisy si- 
multaneous measurement of ji, j 2 , and jz is outlined in 
Fig. 5. The two beam splitters BS provide three copies of 
the original beam by mixing with the vacuum, which are 
directed to three detectors D that are essentially of the 
form of the j\ measuring scheme in Fig. 2. The transfor- 
mations t/2,2 and C^3 )2 placed in front of them transform 
the measurement of ji into measurements of js and j 2 . 
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and Y become the Cartesian coordinates of the plane tan- 
gent to the Poincare sphere at point (j)/(jo) — (0,0, 1) 
|2ll ]. In the same conditions we can consider the bright 
limit of the 12-port detection scheme in Sec. IV, where 
mode a\ plays the role of the local oscillator. When the 
local oscillator is in a coherent state \a) with a — > oo we 
have 
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1+r 2 



of 



(jih) = (jail) 



(AX)' 

(Ay) 2 

= ^{(xy+yj)}. 



l+r z 



(5.7) 



FIG. 5: Illustration of a practical scheme for the noisy si- 
multaneous measurement of ji, j'2, and j'3 in a polarimetric 
context. 



V. BRIGHT LIMIT 

Focusing on the bosonic realization, when the state of 
one of the modes is known, the above measuring schemes 
provide information about the statistical properties of the 
other mode. Let us examine this issue by considering for 
definiteness that the system state factorizes p = p\ ® P2, 
being p\ = \ot)(a\, where \a) is a coherent state a\\a) = 
a\a) with real a for simplicity. In such a case we have 

(h) = a{X), (h) = a{Y), (h) = ±(a 2 -(n)), 

(5.1) 

where X and Y are the quadrature operators of mode a 2 , 
1 



X = 



A 



a 2 



Y = 



A 



a-2 



(5.2) 



and n — a 2 a 2 1S t ne number operator. Concerning the 
covariance matrix we have the following exact series in 
powers of a, valid for any a, 



M = a?M 2 + aMi + M , 



(5.3) 



with 



M, = 



(AX) 
I<(XY + YX)) 




Mi 



(X)(Y) 




1{(XY + YX)) 
(AY) 2 




f(X)\ 
f(Y) 
/ 



(X)(Y) 



with f(A) = (A) + 2(n)(A) - ((An + nA)) and 



Mo = 







(An) 



(5.4) 
(5.5) 

(5.6) 



In the bright limit a — > 00 the a-leading term in 
Eq. (|5.3[) is M 2 so that the angular-momentum covari- 
ance matrix becomes essentially the quadrature covari- 
ance matrix of mode a 2 - This is because in such a limit X 



This agrees well with Eq. (|5.4p . as well as with previous 
analyses of the bright limit of multi-port homodyne de- 
tection [22j . Analyses of homodyne detection with finite 
local oscillator are also available [231 ] . 

Finally we point out that these results are valid in the 
general case beyond two-mode bosonic realizations. This 
holds via the idea of group contraction that applies when 
the state of the system remains in a small enough region 
of the Poincare sphere that can be well approximated by 
the tangent plane [l4T |. 



VI. CONCLUSIONS 

We have provided some general schemes for the practi- 
cal determination of angular-momentum covariance ma- 
trices in different contexts. They allow us to deter- 
mine a global and SU(2)-invariant characterization of 
angular-momentum fluctuations via principal compo- 
nents. This can be of interest, for example, for unam- 
biguous, reference-free characterization of SU(2) squeez- 
ing with applications in quantum metrology, detection of 
many-body and continuous- variable entanglement, and 
light-mediated detection of atomic-spin correlations. 

In particular, we have shown that this task can be ac- 
complished in an exact and simple way by noisy simulta- 
neous measurement of angular-momentum components. 
This provides complete information via the measurement 
of just four observables, instead of the six observables re- 
quired by the general method. Moreover, they are mea- 
sured in a single experimental arrangement. 

It is worth stressing the simplicity of this scheme. The 
minimum number of measured observables required to 
determine the principal variances is three provided we 
know in advance the principal components. Just by 
adding a fourth measurement we no longer need to know 
the principal components, since we can gather enough 
information to determine the whole covariance matrix. 
This includes at once complete information about prin- 
cipal components, principal variances, and the variances 
of angular-momentum components along arbitrary direc- 
tions. 
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